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Abstract. We have studied the statistical properties of 
the energy of individual acoustic modes, extracted from 
310 days of GOLF data near the solar minimum. The 
exponential distribution of the energy of each mode is 
clearly seen. The modes are found to be uncorrelated with 
a ±0.6% accuracy, thus supporting the hypothesis of sto- 
chastic excitation by the solar convection. 
Nevertheless, the same analysis performed on the same 
modes just before the solar maximum, using IPHIR data, 
rejects the hypothesis of no correlation at a 99.3% confi- 
dence level. 

A simple model suggests that 31.3 ± 9.4% of the energy of 
each mode is coherent among the modes studied in IPHIR 
data, correponding to a mean correlation of 10.7 ± 5.9%. 

Key words: Sun: oscillations; magnetic fields - methods: 
analytical, statistical 



1. Introduction 

The quality of GOLF data offers a unique opportunity 
to investigate the excitation mechanism of low degree p 
modes. In the region of 4-6 minutes, several modes can be 
extracted and analysed one by one with a short enough 
time resolution (~ 1.4 days) over one year, allowing an 
accurate statistical treatment of the data obtained. 
Here we deal with the "superficial" energy per unit of 
mass E{t)/M, associated with the mode at the surface of 
the sun. We do not address the issue of the "global" en- 
ergy of the mode, which is related to the superficial energy 
through the shape of the eigenfunction inside the sun. 
Woodard (1984) first revealed the exponential nature of 
the distribution of spectral power in each frequency bin, 
using ACRIM data. 

Goldreich & Keeley (1977) first considered the excitation 
of p modes by turbulent motions near the surface of the 



convection zone. Using the analogy of a damped oscil- 
lator excited stochastically, Kumar, Franklin & Goldre- 
ich (1988) described analytically how the theoretical dis- 
tribution of energy, averaged over a given time- window, 
should depend on the damping time of the mode. Com- 
parisons with real data were performed by Toutain & 
Frohlich (1992) with 160 days of IPHIR data. They found 
in particular that the damping times deduced from the 
linewidths were compatible with those expected in the 
model of stochastic excitation. 

Chang (1996) pointed out that strong localized peaks in 
the energy variations of a damped oscillator excited sto- 
chastically do not necessarily correspond to a strong exci- 
tation, but rather to an exceptional coherent addition of 
random phases. 

This is true provided that the number of independent exci- 
tations per damping time of the mode is large enough, like 
in the case of excitation by solar granules. In this context, 
two different modes excited by exactly the same sources 
would have uncorrelated energies. 

Different modes, however, would be correlated, if the ti- 
mescale between two excitations from a common source 
were longer than their damping time. 
The correlation between the modes energies is therefore 
directly related to the characteristics of their source of ex- 
citation. 

Using 160 days of IPHIR data at the end of 1988 (just 
before the solar maximum, ~ 1990), Baudin et al. (1996) 
concluded that the p modes were likely to be correlated. 
An anticorrelation between the mean solar magnetic field 
and the p-mode power was found by Gavryusev & Gavryu- 
seva (1997) in IPHIR data, while no clear correlation has 
been detected yet in GOLF data (Baudin et al. 1997). 
After checking in Sect. |^ the exponential nature of the dis- 
tribution of energy of p modes in GOLF data, we address 
the issue of their relative independence in Sect. ^, using 
statistical tests based on Montecarlo simulations. 
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These same tests are used to re-analyse IPHIR data in 
Sect, g 

2. Time evolution of the energy of a single mode 

2.1. Method of extraction of the energy 

The energy integrated over a time interval, i.e. the power 
of the mode, was computed by Chaplin et al. (1995) using 
a Fourier transform over short subseries. More sophisti- 
cated methods based on the wavelet analysis were devel- 
opped by Baudin, Gabriel & Gibert (1994) in order to 
analyse the variations of power both with time and fre- 
quency. 

Frequency resolution is not required for our study. Since 
the distribution of energy is likely to be mathematically 
simpler than the distribution of power (Kumar, Franklin 
& Goldreich 1988), we have prefered to extract the energy 
directly. 

Let v{t) be the oscillatory velocity (e.g. integrated over the 
surface of the sun), filtered in the Fourier domain through 
two windows of width Aw centred on the eigenfrequencies 
zbwQ. Its Fourier transform v{uj) is therefore equal to zero 
out of these windows. The time evolution of the energy of 
this isolated mode can be obtained by a bivariate spectral 
analysis, as in Toutain & Frohlich (1992). Here we favour 
a simpler method based on the inverse Fourier transform 
fvit) of the line, translated around w = 0. It is shown in 
Appendix A that the energy of this mode can be written 
as follows: 



fv{t) = 



(1) 
(2) 



This approach is equivalent to the one used by Chang & 
Cough (1995), by means of the Hilbert transform H{v) of 
the velocity, since 4|/„(i)|2 = v'^{t) + H{vY{t). 
If the distribution of velocities is gaussian, then the real 
and imaginary parts of the function fv{t) are two inde- 
pendent gaussian distributions with identical amplitudes 
and variances. Thus Eq. (^ directly implies that the dis- 
tribution of the energy is exponential, as expected. 
Eq. (|l|) shows clearly that the time resolution St of the 
energy, reconstructed by Eq. (|^), is related to the size 
Ai^ = /S.UJ /2tt of the filtering window: 



5t = —. 



(3) 



Denoting by T the total length of the observation, the 
frequency resolution of the Fourier transform is 1/T, and 
the filtering window Av contains p = TAv points. The 
inverse EFT algorithm is used to compute Eq. (nl) and 
define the energy at p successive instants. Eq. (j^ then 
guarantees that the resulting energy is not oversampled. 
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Fig. 1. Time evolution of the energy of the modes 17 < 
n < 25, / = (above) and I = 1 (below). The energy of 
each mode is normalized to its mean value. 



2.2. Application to the GOLF data 

We have considered the set of p modes corresponding to 
17 < n < 25, I = and 1, between 11th April 1996 and 
14th February 1997 (a publication concerning the calibra- 
tion procedure is in preparation). The Fourier transform 
of the resulting velocity over these 310 days allows a fil- 
tering window size of Av = SfiMz {St ~ 1.45 days) for 
this set of modes. The window is symmetric with respect 
to the centroid of the line, wo, which is determined ac- 
cording to Lazrek et al. (1997). The two m-components of 
the mode 1 = 1, however, are not separated. In contrast 
with IPHIR, the width of the window is determined by 
the proximity of another mode (I = 2,3), rather than by 
the level of noise which is here very low. 
Fig. |l| shows the time evolution of the energy of the 18 
selected modes I = and / = 1, normalized to their mean 
energy. The GOLF instrument was stopped for one day 
on 8th September 1996. Four days of signal were removed 
from our statistical study (around the 156th day on Fig. |l]) 
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tio Vp of the estimated variance (denoted by a^) to the 
estimated mean value squared mi: 



Fig. 2. Histogram (20 bins) of the energy of the mode 
I — 0, n = 21, and its cumulative distribution, compared 
to a theoretical exponential distribution. The variance test 
P\v\ compares the observed variance to the theoretical 
one, while the Kolmogorov-Smirnov test Pks depends on 
the maximum distance between the theoretical and the 
observed cumulative distributions. 



in order to account for the stabilisation of the instrument. 
The resulting sample is made up of 210 points. 

2.3. Statistical tests 

Following the picture of a thermodynamic equilibrium be- 
tween the random motions of the convective cells and the 
oscillating cavity (Goldreich & Keeley 1977), we wish to 
compare the observed sample of energies ^i, 1 < i < p, 
with an exponential distribution. Any exponential distri- 
bution is defined by a single parameter, its mean value m. 
Fig. 1^ shows a typical histogram and cumulative distribu- 
tion for the modes extracted from the GOLF data (the 
cumulative distribution is defined as the primitive of the 
density of probability, it increases monotonously from 
to 1). They are compared to an exponential distribution 
whose mean value is estimated from the sample of p 
points. Using the Maximum Likelihood approach, the best 
unbiased estimator of m for an exponential distribution is 
the following: 



(4) 



(i) The variance test 
A simple test consists in checking that the first moments 
of the distribution (mean value and variance) are compat- 
ible with those of a theoretical exponential distribution. 
The variance of an exponential distribution coincides 
with the square of its mean value. We check this property 
by computing, for each mode of the GOLF data, the ra- 



p-i {^i p J2i 



(5) 



Each value is interpreted owing to the cumulative distri- 
bution Py of Vp, obtained if Vp were built from a true 
exponential distribution. Py is computed numerically us- 
ing a Montecarlo method, with 10^ exponential samples 
of p points. For each of the modes selected, Py is the frac- 
tion of these 10^ trials leading to a value of Vp larger than 
the one observed. Since we are interested only in know- 
ing whether the observed Vp is typical of an exponential 
distribution or not, we shall give equal importance to the 
lowest and highest values of the variance by measuring the 
quantity P|y| = 2 min(Py, 100 - Py). 

(ii) The Kolmogorov-Smirnov test 

While the variance test depends only on one partic- 
ular moment of the observed distribution ^i, a more 
global comparison is achieved with the Kolmogorov-Smir- 
nov (KS) test on the cumulative distribution Sp{^). This 
test measures the maximum distance d{Sp,Pm) between 
S'p(^) and a theoretical exponential cumulative distribu- 
tion Pm(x) = 1 — e"^/™. If the mode energies were expo- 
nentially distributed, the statistics of ^(5*^, Pm) would be 
described by a cumulative distribution denoted by P^^- 
Since the mean value m of the reference ditribution is esti- 
mated from the data, we cannot use the standard formulae 
(Numerical Recipies 1992, Chapt. 14.3) to fit P|-g. 
Instead of doing this, we have used a Montecarlo method 
of 10^ samples in order to define the cumulative distribu- 
tion P^g of the distance d{Sp, Pmp). Pks therefore indi- 
cates the fraction of these 10^ trials leading to a distance 
d{Sp, Prrip) larger than the value observed 
For each of the modes selected 
0% would indicate that the observed distribution is too 
far from the theoretical one. A value of P^s close to 100% 
is just as improbable, but would indicate an exceptionnal 
agreement between the theoretical distribution and the 
observed one. 

(iii) Autocorrelation of the artificial exponential distri- 
butions used in the Montecarlo method. 

All of the p modes selected are autocorrelated over a 
timescale comparable to their damping time (2 to 4 days), 
usually deduced from the Full Width at Half Maximum 
(FWHM) of their lorentzian fit in the Fourier space. For 
the sake of accuracy, we have therefore used exponential 
distributions with comparable autocorrelation in order to 
compute the theoretical cumulative distributions Py and 
P^s our Montecarlo simulations. Each one is obtained 
by first creating a time series of a damped oscillator ex- 
cited by a Gaussian noise, and then extracting the energy 
with the method described in Sect. 2T. The damping time 



a value of P^s close to 



of the oscillator is chosen such that it corresponds to a 
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Fig. 3. Result of the Kolmogorov-Smirnov test (square) 
and variance test (plus) of the modes I — and 1 = 1, 
17 < n < 25. As a reference, the horizontal dashed lines 
delimit the upper regions which should contain 68.3%, 
95.4%, 99.7% of the events, respectively. 

FWHM of lfj,Hz in the Fourier space. 
The output of the tests, however, is only slightly modified 
if distributions made of independent points are used. 
For both tests. Fig. ^ shows a very good agreement for 
the set of modes selected. As an exception, the energy of 
the mode n = 23, / = 1 is not exponentially distributed 
{P\v\ = 2.5%, Pks = 1-20%). 

Although the global shape of this distribution can be made 
compatible with an exponential distribution by adopt- 
ing a mean value 10% smaller than the estimated value 
(Pks — 30% for m = 0.9mp), its variance is too large to 
be reconciled with the variance of an exponential distri- 
bution. 

We have also analysed the distribution build with the 18 
modes altogether (each mode is normalized by its mean 
energy). Even with this improved statistics of 18 points, 
the variance and KS tests have not detected any significant 
deviation from an exponential distribution {P\v\ — 60.7%, 
Pks = 72.8%). 

3. Correlation of the individual modes 
3.1. Correlations two by two 

No striking general correlation appears when looking at 
the set of 18 modes displayed on Fig. |l]. Nor does it stand 
out from the computation of the correlations of these 
modes, two by two. Although some large correlations are 
measured (-1-17% between the modes n = 20, ^ = 1 and 
n = 23, 1 = 0), even larger anticorrelations are also found 
(—20% between the modes n = 21, I = 1 and n = 24, 
/ = 1). No general trend is visible, the mean value being 
+0.17% (Fig. I). 

The statistical error of the estimator of the correlation 



Fig. 4. Correlation coefficients of the modes 17 < n < 25 
observed by GOLF, two by two. For each value of n, the 
long and short ticks correspond to Z = and I = 1, re- 
spectively. The symbol + (resp. — ) is used for a positive 
(resp. negative) correlation. The smallest symbols corre- 
spond to correlations smaller than the statistical error, 
intermediate and big symbols correspond to correlations 
smaller than 2 and 3 statistical errors respectively. 



coefficient between exponential distributions, from a 
sample of p points, scales like For our sample 

of 210 points, no effect smaller than 7% can therefore 
be detected. Altogether, 5.5% of the couples (17 out of 
306 couples) present correlations contained, in absolute 
value, between 2 and 3 standard deviations, which is not 
very significant (4.3% would be expected for a normal 
distribution of the statistical error) . 

Nevertheless, a more sensitive indicator can be construc- 
ted, in order to determine the mean correlation coefficient 
more accurately. 



3.2. Test of the null hypothesis. Comparison with a 
Gamma distribution 

If k distributions are independent (null hypothesis), the 
variance of their sum should be equal to the sum of their 
variances. This test was performed by Baudin et al. (1996) 
with IPHIR data, who normalized the distribution by their 
level of noise, and found some discrepancy. 
A fundamental feature of our method is the use of the ex- 
ponential nature of each individual energy distribution, in 
order to compute the standard deviation of our estimate 
of the correlation, and therefore the confidence level of our 
conclusions. 

We denote by Tk,p the sum of k distributions of energy, 
made of p events, where each of the distributions is nor- 
malized by its estimated mean energy. Since each distribu- 
tion appears to be exponential within the statistical error 
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(Section 2), T^.p ought to resemble a Gamma-distribution 
of order k (denoted by r^) if they are independent, or an 
exponential distribution of mean value k if they are all 
identical. The null hypothesis can therefore be tested by 
comparing the observed distribution Tk,p with the theo- 
retical Ffc-distribution, using the variance and KS tests. 
Denoting by Cij the correlation coefRcient between the 
modes i and j, and C their mean value, the variance of 
Tfc = hmp^oo Tfc^p is directly related to these correlations: 



var(Tfc) = fc + 2 ^ C,,j-. 

l<i<i<fc 

= k + k{k- 1)C. 



(6) 
(7) 



If the modes are independent, the standard deviation cr^j^ 
of the variance estimator of ~ Tk is 



k 

T" 

p2 



(8) 



Consequently, another way of testing the null hypothesis 
is the comparison of the variance of T^.p with yar(Tk) = fc, 
in units of the statistical error ir^ar- 

Even if the k distributions defining Tk.p are independent 
and exponential, an additional error of the order of p^^/^ 
is introduced in the estimation of the mean value of each 
exponential distribution. We use a Montecarlo method of 
10^ trials made from independent exponential distribu- 
tions, in order to define the cumulative distributions Py^ 
for the outcome of the variance test, and P^'g for the 
outcome d{Tk,Tk,p) of the KS test. Here again, we have 
used autocorrelated exponential distributions in the Mon- 
tecarlo simulations. 

Eq. (|^) indicates that the variance of the distribution 
gives a direct measure of the mean correlation among the 
modes: 



C = 



var(Tfc^p) k 
k{k~l) 



O 



k{k - 1) 



(9) 



This formulae, however, is not directly useful without an 
expression of the statistical error tTvar associated with the 
estimator of the variance. If the modes are correlated, 
computing it requires some additionnal ass um ptions about 
the properties of the correlation (Sect. ^^). Neverthe- 
less, (Tvar coincides with Eq. (H) to first order. Together 
with Eq. (H), the smallest correlation detectable with this 
method scales as follows: 



C ■ = 



fc(fc- 1) 



(10) 



which is a factor -^2 / k smaller than the sensitivity of the 
correlation coefficient two by two. A better sensitivity is 
therefore obtained by summing a large number of modes. 
However, the hypothesis of a constant correlation between 
the modes might be questionable if the range of frequen- 
cies is large, especially since the mean energy and the life- 
time of the modes vary significantly with frequency. 




Time (days) 



Fig. 5. Time evolution of Tg ^io, the sum of the normal- 
ized energies of the 9 modes ^ = 0, 17 < rt < 25, observed 
by GOLF. 




Fig. 6. Histogram (20 bins) of the energy and cumulative 
distribution of the sum of the 9 modes I — 0, 17 < n < 25 
observed by GOLF, compared to a Fg distribution. 



3.3. GOLF results 

Both tests are of course very sensitive to the presence of 
a gap in the data. If x % of the data were filled with zeros 
due to an interruption of the instrument, the variance of 
'^k,p would be increased by a factor ~ (1 -|- k)x %. Conse- 
quently, we have carefully removed from our samples the 
points corresponding to these gaps. 

The sum of the energies of the 9 modes I — 0, 17 < n < 25, 
normalized to their mean energy, is shown in Fig. ^j. We 
note in passing that the clear gap in the data appearing 
around the 156th day confirms the validity of our proce- 
dure for the extraction of the energy. As before, four days 
of signal have been removed from our statistical study to 
account for the stabilisation of the instrument, resulting 
in a sample made up of 210 points. Their distribution 



6 



T. Foglizzo, R.A. Garcia et al. : Are solar acoustic modes correlated ? 



is successfully compared to a Fg-distribution in Fig. || 
(P|y| 49.9% and Pks = 76.2%). The same test per- 
formed on the 9 modes ^ = 1, 17 < ti < 25, obtains 
P\v\ = 31.6% and Pks = 34.2%. Applied to these 18 
modes altogether, the tests confirm again the null hypoth- 
esis {P\v\ = 95.8% and Pks = 66.5%). 
The correlation being low, we may use Eq. (|^) with the 
statistical error given by Eq. (^), and obtain a confir- 
mation of the absence of correlation, with an error bar: 
C = 0.7 ± 1.4% for 9 modes Z = 0, and C = -0.1 ± 0.6% 
for 18 modes I = 0,1. 



3.4- Test of the "X-hypothesis". Correlation due to an 
additive common signal 

A refined estimate of the correlation can be obtained by 
making some assumptions about its origin. We build in 
Appendix B a simple model where the excitation is a 
mixture of two types of sources, which, taken separately, 
would result in uncorrelated/highly correlated modes en- 
ergies respectively. 

The first type represents the granules, which produce so 
many excitations per damping time that the correlation 
among the modes energies is close to zero. 
The second type is hypothetical. It could be produced by 
some isolated events, possibly of magnetic origin, sepa- 
rated by a time comparable to or larger than the damping 
time of the modes considered. 

We assume that the mode response to an excitation is lin- 
ear, and therefore the response to a mixture of sources is a 
superposition of the answers to the two types separately. 
Our model depends on a single parameter A, namely the 
fraction of the energy of each mode due to the second type 
of sources. 

We define in Appendix B the theoretical distribution func- 
tion F^ corresponding to such correlated modes energies. 
If the model is applicable, the distribution Tk,p ought to 
converge, for p oo, towards a well defined distribution 



denoted by F^, such that: 
Tt°ix) = Ffc(x), 



F 



-(.) . -e- 



(11) 
(12) 



The variance and KS test can therefore be used, for vari- 
ous values of A, in order to test this "A-hypothesis" . 
With the definition of Appendix B, the correlation coeffi- 
cient C is related to the coefficient A as follows: 



C = A" 



(13) 



The statistical error dvar associated with the estimator of 
the variance also depends on A according to Eq. ( |B.6D . 
Eq. (|^) can then be used to determine the mean correla- 
tion C, with a consistent statistical error. 
We do not expect higher order moments of the distribu- 
tion Sk to be more sensitive to a correlation between the 
modes, since we prove in Appendix B that they also vary 




c (%) 



Fig. 7. Comparison of GOLF data {p — 210 points) with 
a theoretical distribution made of k modes with a uniform 
correlation C, using the variance and KS tests. The vertical 
dotted line delimits the sensitivity limit {Cmin — 0.6%) 
defined by Eq. ( |lO| ) for fc = 18 modes. It coincides with 
the value below which the corresponding variance and KS 
tests remain inside the upper 68.3% region. 



like A^ at first order. 

We also demonstrate that the shape of the cumulative dis- 
tribution varies like A^. The sensitivity limit of the KS test 
is therefore expected to be comparable to the sensitivity 
limit of the variance test. 

Here again, normalizing the distributions by their esti- 
mated mean value introduces a bias, which we take into 
account using a Montecarlo method. For each value of 
k, A considered, we compute from 10^ trials the theoreti- 
cal distribution F^, and use 10^ other trials to define the 



cumulative distributions p^''^^^ 



and Pks'^ which are used 



for our variance and KS tests. For the sake of simplicity, 
the effect of the autocorrelation of each mode is neglected 
here. Indeed, we know from Sect. 2.3 and 3.2 that it in- 



troduces very small corrections on the cumulative distri- 
butions Py^ and Pks- 

We define the error bar of the correlation coefficient as the 
range of values of C = A^ within which the test {P\v\ or 
Pks) remains inside the upper 68.3% region, by analogy 
with the statistics of normal distributions. 
Fig. 1^ shows that the variance and KS tests, applied to 
18 modes for various values of A, stays within the upper 
68.3% region for C < 0.6%, which coincides with the sta- 
tistical limitation expressed by Eq. ([lo|). 
Nevertheless, while the measurements made by GOLF are 
compatible with a total lack of correlation of the modes, 
our tests cannot exclude that up to A„iax ~ 8% of the 
energy is common to the modes. 



4. Comparison with IPHIR data 



p,,(n) P|viC+) 

+ □ 
□ □ 
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Frequency (/aHz) 



Fig. 8. Variance test (plus) and KS test (square) for the 
11 modes / = 0, 19 < n < 23, and ^ = 1, 18 < n < 23 
observed by IPHIR with a fihering window of 6 /iHz. 



4.I. Extraction of 11 modes 

Since the conclusions of Baudin et al. (1996) about 160 
days of IPHIR data were obtained from time variations 
of the power instead of the energy, using a different nor- 
malization and a different time resolution, we have first 
re-analysed these data with the method described above, 
using the same 11 modes (^ = 0, 19 < n < 23, and ^ = 1, 
18 < n < 23). The central frequency ojq is taken from 
Toutain & Frohlich (1992). The higher level of noise lim- 
its the size of the filtering window to 6/iIIz, leading to a 
time resolution of 1.9 days, and a statistical study on 82 
points. We have removed the data surrounding two gaps 
in the series, around the 5th and the 61st day. The result- 
ing sample is shortened to 78 points only. 
According to Fig. ||, the distribution of energy of each of 
the 11 modes of IPHIR is compatible with an exponential 
distribution (apart from the mode I = 1, n = 22). 
The correlation coefficient of the modes energy, two by 
two, is shown in Fig. ^. The mean value is 3.9%, the statis- 
tical error being 11.3%. Altogether, 4.4% of the couples (2 
out of 45 couples) present correlations contained, in abso- 
lute value, between 2 and 3 standard deviations, which is 
comparable to the 4.3% expected for a normal distribution 
of the statistical error. 

4-2. Test of the null hypothesis 

We have also extracted the same modes, with the same 
filtering window, from the first 153 days of GOLF data to 
obtain a comparable sample of 78 points. The difference 
between the two series appears on the distribution Tn^yg 
shown in Fig. |l^, where both the variance and the KS tests 
indicate that the modes are likely to be less independent 
in IPHIR data than in GOLF data. 

The tests applied to GOLF data are compatible with the 



Fig. 9. Correlation coefficients of the 11 modes observed 
by IPHIR, two by two. For each value of n, the long and 
short ticks correspond to ^ = and ^ = 1, respectively. The 
symbol + (resp. — ) is used for a positive (resp. negative) 
correlation. The smallest symbols correspond to correla- 
tions smaller than the statistical error, intermediate and 
big symbols correspond to correlations smaller than 2 and 
3 statistical errors respectively. 



null hypothesis (P|v| = 35.0%, Pks = 82.8%), which is 
consistent with the results of Sect. |3[ 
By contrast, the same tests applied to IPHIR data reject 
the null hypothesis with a 99.32% confidence level with 
the variance test, and a 95.72% confidence level for the 
KS test (P|y| = 0.68%, Pks = 4.28%). 

4-. 3. Test of the "X-hypothesis" 

While the cumulative distribution Tnjs of GOLF does 
not show any systematical trend when compared to Fn, 
the cumulative distribution Tn^rg of IPHIR shows a clear 
trend. This trend is successfully suppressed when com- 
pared to the distribution T^i^ (Fig. |l^). Fig. |ll| shows 
that within our simple model, the signal of IPHIR would 
be absolutely normal as regards our tests {P\v\ = 85.0%, 
Pks = 42.6%) if a fraction A = 25% of each mode energy 
were common to all the modes, correponding to a mean 
correlation C = 6%. 

Error bars are obtained by varying the parameter A: the 
variance test leads to C = 6.1 ± 3.3%, while the KS test 
obtains C = 10.7 ±5.9%. 

Moreover, the correlation computed from Eq. (m with the 



statistical error given by Eq. (B.6) is C = 5.0 ± 3.6%. Al- 
though this analytical estimate is less reliable than the 



tests based on Montecarlo simulations (Eq. (B.6) neglects 



the error in estimating the mean energy of each mode) , it 
is useful as a quick check of the results. 
It is therefore comforting to notice, as can be seen in 
Fig. 11, that the range of correlations defined by these 
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Fig. 10. Histogram (15 bins) of Tn.yg and cumulative 
distribution, compared to a Fn distribution for GOLF 
(above) and IPHIR (below). In the IPHIR case, the dashed 
lines correpond to the theoretical distribution F^j^-^^. 

three methods overlap in the range 4.8% < C < 8.6%, 
correponding to 21.9% < A < 29.3%. Of course, this over- 
lapping region cannot be directly interpreted in terms of a 
standard deviation. We shall adopt the conservative range 
obtained with the KS test, which takes the full distribu- 
tion into account: C ~ 10.7 ± 5.9%, corresponding to a 
fraction A = 31.3 ± 9.4%. 

4- 4- Additionnal checks 

In order to check the possibility that the correlation might 
come from a multiplicative noise (such as due to a point- 
ing noise), we have computed the correlation between 11 
windows of noise centered 20fiHz (resp. 21 ^Hz) to the 
right of each mode. This test indicates that the noise it- 
self is not correlated, with P\y\ = 77.2% and Pks = 25.5% 
(resp. P\v\ = 68.4% and Pks = 34.3%). 
We have checked the effect of changing the size of the fil- 
tering window to 4 /iHz (no noise, but low statistics of 
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Fig. 11. Comparison of IPHIR data [p — 78 points) with a 
theoretical distribution made of fc = 11 modes containing 
a fraction A = C^/"^ of energy in common. The vertical 
dotted lines delimit the range of the correlation coefficient 
determined by Eq. (^. The KS test remain inside the 
upper 68.3% region for C — 10.7 ± 5.9%, corresponding to 
a fraction A = 31.3 ±9.4%. 

52 points) and 8 /LtHz (good statistics of 106 points, but 
IPHIR is influenced by the noise). While a smaller filter- 
ing window still favours A ~ 25%, a larger window takes 
into account a significant fraction of uncorrelated noise, 
as expected, resulting in a slightly lower value of A. 
One might also suspect that the discrepency between the 
IPHIR distribution Tn 73 and a Fn distribution is due to 
the mode Z = 1, n = 22 which is not well fitted by an expo- 
nential distribution (see Fig. Nevertheless, performing 
the same analysis without this particular mode leads to 
the same conclusion: P\v\ = 0.28% and Pks = 6.9% if 
A = 0, while P\v\ = 78.2% and Pks = 38.5% if A = 25%. 

5. Conclusion 

The exponential nature of the energy distribution of each 
mode has been used to compute their mean correlation 
coefficient with a consistent error bar. 
Two tests based on Montecarlo simulations, and one ana- 
lytical formulae, applied to 310 days of GOLF data, sup- 
port the null hypothesis of no mean correlation among 
the modes 17 < n < 25, ^ = 0, 1, with an accuracy of 
C = ± 0.6%. 

Our analysis of the modes correlation in IPHIR data, using 
these statistical tools, gives an accurate statistical support 
to the tentative conclusions of Baudin et al. (1996). The 
variance test rejects the null hypothesis with a 99.3% con- 
fidence level. 

The presence of a clear correlation among p-mode energies 
in IPHIR data strongly constrains the standard picture of 
stochastic excitation. If really of solar origin, it suggests 
the existence of an additional source of excitation, other 
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than the granules. We have built a one parameter model 
of random excitations separated by a time comparable to 
the damping time of the modes, added to the usual gran- 
ule excitations. IPHIR data are fully compatible with this 
"A-hypothesis" , if a fraction A = 31.3 ±9.4% of each mode 
energy is due to this additionnal source of excitation, re- 
sulting in a mean correlation C = 10.7 ± 5.9% among the 
modes. 

On the other hand, the absence of correlation in GOLF 
data support the standard picture of stochastic excitation 
by the granules only. 

This difference between IPHIR and GOLF data can be 
interpreted as a change from A ~ 31.3 ± 9.4% in IPHIR 
data to less than 8% in GOLF data. 

This evolution could be related to the change in magnetic 
activity, since the GOLF data correspond to a period close 
to the solar minimum while the IPHIR data correspond 
to a period closer to the solar maximum. If this is true, 
a confirmation will be obtained by performing this same 
analysis on GOLF data when we approach the solar maxi- 
mum, in a couple of years. VIRGO data will also be useful 
in order to identify the possible role of the measurement 
techniques (velocity/intensity) in the determination of the 
correlation. 

However, the mechanism by which the magnetic field in- 
fluences the excitation of the modes, i.e. the nature of 
these hypothetical exciting events remains to be explored 
in more detail. 
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Appendix A: Extraction of the time evolution of 
the energy 

Let us consider the real displacement y{t) filtered through a 
double window with a width Alj, centred on the frequencies 
±ijJo: 



y{t) 



y{ijj)e'^ Auj + c.c. 



(A.l) 



where c.c. is the complex conjugate of the first term. By anal- 
ogy with an oscillator of eigenfrequency ujo, the energy is de- 
fined as the sum of a kinetic and a potential part; 



/ ,\ _ 1 / 2 , 2 2x 



(A.2) 



From Eq. (|l|) and (A.l), the filtered displacement and velocity 
can be written as follows: 



y{t) = 7^eal{2e'"''Vy(^)}, 
v{t) = TZe&\{2 e'"°*/4^)} , 



(A.3) 
(A.4) 



E_ 

M 



it) = i/„r+a;o'i/«r 

~f7^eal{e^'"«*(/„^(^)+^c'^'(^))} 



We first deduce from the relation v = dy/dt that and fy are 
related as follows: 



fv{t) = icjo / y{uJo + (^)e'^* + —'] duj, 



^.ofyit){l + o{^)} 



(A.6) 



(A.7) 



From Eq. (|l|), the non-zero Fourier components of fy (and fv) 
correspond to frequencies between —Auj/2 and +Auj/2. If, as 
is usually the case, Au <^ uuo, the high frequency oscillations 
of exp it2iLoot are well se parat ed from the slower variations of 
fy{t) and fv{t), and Eq. (A..5) is approximated by Eq. (bl). 



Appendix B: One parameter model of correlated 
modes 

The amount of energy which is coherent among the modes can 
be estimated by constructing a simple one-parameter model as 
follows. Indexing the modes by j £ {1, fc}, we assume that the 
velocity residual of each mode Vj{t) = Vj{t) + ajVo{t) is made 
of a superposition of two independent signals, where vo{t) is 
common to all the modes, and all the Vj{t), j £ {0, fc}, are 
independent. 

Using the filtering method described in Section 1, we introduce 
the parameter 9j as: 



/v,(t) = fv,it) + a,,U„it) 



+ sin ( 



(B.l) 
(B.2) 



where {rj,ij), j £ {0, fc}, are independent normalized normal 
distributions. The quantity Xj = sin 6j can be interpreted as 
the ratio of the energy in the common signal to the total energy 
of the mode. We denote by ej the energy of the signal filtered 
in the Fourier space, normalized to its mean value: 

Cj = - ^(rj cosOj + rosinOj)^ + {ij cos9j + iosinOj)'^^ . (B.3) 



(B.4) 



For the sake of simplicity, Xj = A is assumed to be indepen- 
dent of the mode j. This is equivalent to assuming that the 
correlation is uniform among the modes. 
The sum of the normalized energies is defined as: 



Tfc 



(B.5) 



The variance of the estimator of the variance depends on the 
fourth moment of the distribution, and is equal to: 



2k{k + 3) 



2k[k - 1)A^ 



[2(fc + 9) + 12{k - 2)A + (4fc^ ~ lOfc -I- 9)A^]. (B.6) 



Let us show that the higher order moments of the distribu- 
tion Tfc vary like A^, to first order. 



(A.5) /i, = ((Tfe-fc)'), 



(B.7) 
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where ( ) denotes the expectation value of the distribution. 
As the transformation —9 does not change the distribu- 



tion defined by Eq. (B.3), only the even powers of sin^ can 
contribute to fii. It can therefore be expande d into powers of 
A = sin^ 6. Let us develop the product in Eq. (B.7) and prove 
that the term of order A is zero. We use below the special 
relation between the centred moments jli,jli-\,jli-2 of a Tk- 
distribution of order k: 



fii = fii+AX + OiX"), 




j=i 



l{l ~ 1) {fil-l + fc/ii-2) 

0. 




Let us now show that the probability density (x) of Tfe also 
varies like A'^, to firs t order. 

Given the Eqs. (B.3)-(B.J:) defining T^, we can expand r^(a;) 
in powers of A. 



The moment of order / is defined as; 



{x — fc)'rfc(a;)da; 



fii+\ ix~kyg{x)dx + 04X^). 



(B.12) 

(B.13) 
(B.14) 



The fact that every moment varies at least like A (Eq. B.ll) 
implies: 



{x — k) g{x)dx — 0. 



(B.15) 



The only function g{x) satisfying Eq. ( B.15 ) for any value of 
I is g = 0. Therefore both r^(a;) and its primitive (i.e. the 
cumulative distribution of Tk involved in the KS test) vary 
like A^ to first order. 

We conclude that the sensitivity of the KS test is comparable 
to the sensitivity of the variance test. 
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